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We investigate a spin reorientation transition (SRT) in ultrathin magnetic films by Monte-Carlo simula- 
tions. We assume that the lateral size of the film is relatively small and it has a single-domain structure. To 
gain insights into the SRT, we measure a free-energy as a function of perpendicular and in-plane magneti- 
zations. As a result, we find that the system is in a superparamagnetic state at the SRT temperature. The 
disappearance of magnetization around the SRT temperature, which is observed in experiments, emerges 
due to dynamical fluctuations in magnetization which are inherent in the superparamagnetic state. This 
observation is in contrast to that in large ultrathin magnetic films that the disappearance of magnetization 
is caused by a static magnetic structure with many complex domains. 
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1. Introduction 

In the last several decades, ultrathin magnetic films have 
studied extensively due to both the fundamental interest of 
low-dimensional magnetism and the numerous perspectives 
of applications. These studies have revealed that complex 
magnetic order and curious phenomena are induced in ul- 
trathin magnetic films by competition among several interac- 
tions such as ferromagnetic exchange interactions, magnetic 
anisotropy energies, and magnetic dipolar interactions. A spin 
reorientation transition '~ 5) (SRT) is one of such phenomena. 
This is a transition from a in-plane magnetized state to a per- 
pendicularly magnetized state. This transition is induced by 
decreasing the thickness of the film or by decreasing the tem- 
perature. The SRT has been extensively investigated both ex- 
perimentally and theoretically until now (see reviews 6, 7 and 
references therein). 

The main issue we address in the present study is the 
disappearance of magnetization around the SRT temperature 
Tsrt- 1 ' 4 ' This abrupt drop in magnetization around Tsrt is 
called pseudogap.^ Papas et al. proposed that there are two 
possible origins which cause this pseudogap. The first is a 
dynamical origin that the system is in a paramagnetic state 
because perpendicular and in-plane anisotropics compensate 
with each other around Tsrt. As a result, magnetization heav- 
ily fluctuates with time, and the time average of each com- 
ponent of magnetization, which is measured in experiments, 
becomes zero. The second is a static origin that complex mag- 
netic domains, which cause an almost complete loss in the 
total magnetization of the sample, emerge in the vicinity of 
Tsrt- For large ultrathin films with the lateral size of several 
tens of micrometers, it was observed that a complex domain 
structure appears near the SRT temperature. 3 - 5) Such behavior 
was also observed in Monte-Carlo (MC) study 8 ' if the corre- 
sponding system size is large enough. 9 ' These results strongly 
indicate that the pseudogap is caused by the static origin for 
large films. However, it has not been clarified whether the 
pseudogap emerges due to the dynamical origin even if the 
lateral size of the film is not large and the system has a single- 
domain structure. 

To address this issue, we examine the SRT of ultrathin 
magnetic films with focusing on the pseudogap around Tsrt- 



The corresponding size of the film is chosen to be relatively 
small so that the system has a single-domain structure. To gain 
insights into the properties of the SRT, we measure a free- 
energy as a function of the two order parameters of the SRT, 

1. e., perpendicular magnetization m x and in-plane magnetiza- 
tion ni\\. To our knowledge, this is the first time that such free- 
energy is measured in the study of the SRT. The free-energy 
is calculated numerically by using the method proposed in 
ref. 10. As a result, we found that the free-energy around Tsrt 
has a structure that low free-energy region widely spreads in 
the (m ± ,m\\) space. Because the amplitude of the magnetiza- 
tion m = (nv\ + m 2 ) 1 / 2 is non-zero in the low free-energy 
region, spins are aligned in the same direction. However, the 
direction of the magnetization heavily fluctuates with time as 
a result of the exploration inside the low free-energy region. 
In short, the system is in a superparamagnetic state at Tsrt- 
Our results indicate that the pseudogap emerges due to the 
dynamical origin even in relatively small ultrathin magnetic 
films with a single-domain structure. 

2. Model 

In the present work, we consider a spin model on a 32 X 
32 x 1 square lattice with open boundaries. The Hamiltonian 
of the model is described as 



(i,f> 
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S r Sj (S r r u )(Sj-r u )) 



(1) 



where S, is a classical Heisenberg spin with an absolute value 



of unity, r, 



= r. 



rj, r, is the position vector of a site 



i in the unit of the lattice constant, and ry = |fy|. In the 
right hand side of eq. (1), the first term, the second term, 
and the third term denote ferromagnetic exchange interac- 
tions, dipolar interactions, and magnetic anisotropy energies, 
respectively. The sum of exchange interactions runs over all 
of the nearest-neighboring pairs. We assume that magnetic 
anisotropy is uniaxial and its easy axis is parallel to the z- 
axis, where the z-axis is perpendicular to the film. This type 
of model has been often used in the study of the SRT. 8 ' n ~ 17) 
The three parameters /, Cd, and C u in eq. (1) represent the 
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Table I. The four sets of parameters used in the present study. The mesh 
size <5 is calculated from eq. (2) and the ratio I J by assuming that A = 
30.2 X irr 12 J/m and / s = 1.79 T. The size is calculated from the mesh 
size i5 and the lattice size 32 X 32 X 1 . 





CJJ 


CJJ 


6 [ran] 


size [run 3 ] 


case 1 


0.007 


0.05 


1.44 


46.1 X46.1 x 1.44 


case 2 


0.028 


0.2 


2.89 


92.5 X 92.5 X 2.89 


case 3 


0.056 


0.4 


4.08 


131 x 131 x 4.08 


case 4 


0.070 


0.5 


4.56 


146 x 146x4.56 



strength of each energy term. In the present work, we per- 
formed simulations for four different sets of the parameters. 
The values of Cd/J and C u /J in the four cases are shown in 
Table I. The ratio Cd/C u , which is 0.14 in all of the four cases, 
is chosen so that the SRT is clearly observed. In micromag- 
netic calculations, the three parameters in eq. (1) are calcu- 
lated from the material parameters and the mesh size as 

J = 25 A, C d = J^6 3 /(4-nfi Q ), C u = K a 6 3 , (2) 

where 5 is the mesh size, A is the exchange stiffness con- 
stant in continuous limit, J s is the saturation polarization, K u 
is the magnetic anisotropy constant, and jUo is the vacuum 
permeability. The mesh size 5 shown in Table I is calculated 
from the ratio Cd/J by assuming that A = 30.2 x 10~ 12 J/m 
and J s = 1.79 T. These values are the same as the mate- 
rial parameters of cobalt, 18 ' which is one of the representa- 
tive magnetic materials. The sizes of the film shown in Ta- 
ble I are calculated from the mesh size and the lattice size. 
The magnetic anisotropy constant K u is calculated from the 
ratio C d /C u = 0.14 as 1.45 x 10 6 J/m 3 . Since the magnetic 
anisotropy constant of cobalt is 4.53 x 10 5 J/m 3 , 18) this value 
is about three times larger than that of cobalt. In the present 
study, we examine the mesh-size dependence by comparing 
the data of these four cases. 

3. Method 

We next explain how we measure free-energy as a function 
of m_L and my. The free-energy is defined by 

exp [-{3F 06; m ± , m,,)] = Tr {Sj) exp [-/W{S,}] 

X6 (m_L - mliSj}) 5 (in l{ - rajjfS,}) , (3) 

where is the inverse temperature and < H{Si] is the Hamil- 
tonian defined by eq. (1). The trace in the right-hand side of 
eq. (3) runs over all of possible configurations. Note that the 
right-hand side of the equation is proportional to the proba- 
bility that spin configurations with magnetizations (mj_,W||) 
are observed. »z^{S,} and mjJ{S,} are the perpendicular and in- 
plane magnetizations calculated from a configuration {S,}, re- 
spectively, i.e., 

<{*)=£E,4 (4) 

where N is the number of spins. 

We used the method proposed in ref. 10 to calculate 
F (j3;m ± ,m^) numerically. In this method, a variant of the 
Wang-Landau method 19 ' 20) is combined with the stochas- 
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Fig. 1. (Color online) The temperature dependences of the perpendicu- 
lar magnetization m± (open circles) and in-plane magnetization ma (open 
squares) in the case 1 (upper left), 2 (upper right), and 3 (lower left) in Ta- 
ble I. In the measurement, the temperature was gradually cooled from 2.0 J 
to 0.025 J in steps of AT = 0.025 J. The system was kept for 100, 000 MC 
steps at each temperature. The data of the latter 50, 000 MC steps were 
used to calculate the thermal averages of the two magnetizations. The av- 
erage was also taken over 10 different runs with different initial states and 
random number sequences. 

tic cutoff (SCO) method, 211 which was recently invented for 
long-range interacting systems. This method enables us to 
calculate F (J3;m ± ,m\\) with reasonable computational time. 
Furthermore, long-range dipolar interactions are taken into 
account without any approximations by the use of the SCO 
method. The detailed conditions of free-energy measurements 
are as follows. We set the initial value of the modification 
constant AF in ref. 10, which is related with the modification 
factor / in the Wang-Landau method 19 ' 20 ' by / = exp(AF), 
to unity. We stopped our simulation after we halved AF 20 
times. Therefore, the final AF is 2 -20 . The magnetization 
space (m ± ,m||) was divided into 10,000 bins with the area 
of 0.01 x 0.01, and free-energy was calculated for each bin 
which satisfies the inequalities < (to 2 + to 2 ) 1 / 2 < 0.85. 
The histogram of two magnetizations //(otj_, m\\) was checked 
every 10,000 MC steps. We regarded the histogram as flat 
when H(m ± , m\\)'s for all the magnetizations are not less than 
80% of the average value of the histogram. The SCO method 
was applied only to dipolar interactions. A potential switch- 
ing procedure in the SCO method was performed every 10 
MC steps. We performed the free-energy measurement for 10 
different runs with different initial states and random num- 
ber sequences to estimate the means and errorbars of the free- 
energy. The estimated errorbars were rather small. They were 
smaller than 0.03 J for all of the data we will show hereafter. 

4. Results 

To check whether this model exhibits the SRT or not, we 
first measured the perpendicular and in-plane magnetizations 
as a function of the temperature in the three cases 1, 2, and 
3 in Table I. Figure 1 shows the result. We hereafter set the 
Boltzmann constant k% to unity and use J as a unit of tem- 
perature. In this measurement, the temperature was gradually 
cooled from 2.0 / to 0.025 J in steps of Ar = 0.025 J, and the 
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Fig. 2. (Color online) Free-energy difference F^g in the case 2 is plotted 
as a function of (m ± , mir) for two different temperatures. F^g is defined by 
eq. (6). The left and right panels show the data for T = 7/srt - 0.04 J and 
those for T = 7/srt + 0.04 J, respectively. 7/srt is 0.33 /. In both figures, 
Fdiff is plotted only for a region in which F^g is smaller than 2.5 /. The 
low free-energy region defined by the inequality Fag < T is denoted by a 
solid white line. The average was taken over 10 different runs. 

thermal averages of the two magnetizations were measured at 
each temperature. As the temperature decreases, m\\ starts to 
increase rapidly around T ~ 0.8 J. On the contrary, such rapid 
increase does not occur in m ± around this temperature. There- 
fore, my is dominant at intermediate temperatures. However, 
my starts to drop around 0.4 J, and m ± overtakes my around 
0.3 J. These data clearly show that the present model exhibits 
the SRT We also see that the SRT temperature Tsrt, which is 
defined by the temperature where m ± is equal to ni\\, is around 
0.3 J in all of the three cases. This means that Tsrt mainly de- 
pends on the ratio Cd/C u . Recall that the ratio is fixed to 0.14 
in all the four cases (see Table I). We also performed MC sim- 
ulation in the case 4, and found that a two-domain structure 
emerges in this case. Because the purpose of the present study 
is to investigate the SRT in single-domain magnetic films, we 
hereafter show the results of the three cases 1-3. Since the size 
of the case 3 is close to that of the case 4, the three cases 1-3 
cover almost the whole range of the size in which a single- 
domain structure is retained. 

One may consider that a pseudogap does not exist in this 
model because neither m ± nor my is zero at Tsrt- However, 
this result is not incompatible with the presence of a pseu- 
dogap due to the dynamical origin. It should be recalled that 
each component of magnetization m^ (fi — x,y,z), which can 
be both positive and negative, disappears at T^kt- When spins 
are aligned in the same direction and the direction changes 
with time, the time average of becomes zero. However, 
the time averages of m ± and my do not become zero because 
they are always positive by definition. 

We next measured the free-energy F(J3; m ± , my) defined by 
eq. (3). In Fig. 2, free-energy difference /\jiff in the case 2 is 
plotted as a function of (m ± , nin), where F^g is defined by 

^diffOS; m ± , my) = F(/3; m x , m\\) - F miB (J3), (6) 

and 

f min 08) = min F(j3;m ± , mil ). (7) 

Figures 2(a) and 2(b) show the data for T < Tsrt and those 
for T > 7srt> respectively. We see that F^ for T < Tsrt has 
a global minimum around (m x ,m||) ~ (0.75,0.15). In con- 
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Fig. 3. (Color online) Free-energy difference F^g at Tsrt is plotted as a 
function of (m±, mil) for the case 1 (upper left), 2 (upper right), and 3 (lower 
left). Fdjjj is plotted only for a region in which F^g is smaller than 2.5 J. 
The low free-energy region defined by the inequality F^g < T is denoted 
by a solid white line. The average was taken over 10 different runs. 



trast, when T > Tsrt, ^difr has a global minimum around 
(m ± , my) * (0, 0.6). The location of the global minimum shifts 
discontinuously as the temperature is changed across T^rj- In 
this sense, we can regard the SRT as a first-order transition. 
This discontinuous shift of the global minimum explains the 
reason why m ± and m\\ in Fig. 1 change abruptly around 7srt. 
This result is consistent with the previous results that single 
monolayer magnetic films with a single-domain structure ex- 
hibit a first-order transition. 12 ' 17,22) A similar discontinuous 
shift in the global minimum was observed in the other two 
cases 1 and 3. 

We next focus our attention to the free-energy structure at 
Tsrt. In Fig. 3, we show F&ft for all of the three cases. There 
are several differences in the free-energy structure among the 
three cases. When Ca and C u are small (case 1), the free- 
energy has a circular structure. In the case 2, this free-energy 
structure is slightly deformed into a elliptical one. For large 
Cd and C u (case 3), the deformation proceeds further, and the 
free-energy has two local minima and a free-energy barrier in 
between. However, a common feature among the three free- 
energy structures is that a low free-energy region defined by 
the inequality F^ff < T, which is denoted by a solid white line 
in Fig. 3, widely spreads in the (m ± , my) space. In the cases 1 
and 2, low free-energy region involves both the perpendicu- 
larly magnetized state with my ~ and the in-plane magne- 
tized state with m ± « 0. There is no (or vanishingly small) 
barrier between the two states. In the case 3, there is a free- 
energy barrier like an ordinary first-order phase transition. 
However, the barrier is rather small. Because single-domain 
magnets are known to exhibit large fluctuations in magnetiza- 
tion if the energy barrier is small in comparison with k^T, 2 ^ 
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Fig. 4. (Color online) The amplitude of magnetization m = (m 2 ± + mjj) 1 ' 2 
and the z-component of magnetization m z measured at 7srt are plotted as 
a function of time for the two cases 1 (upper panel) and 3 (lower panel). 



Fig. 5. (Color online) In-plane components of magnetizations m x and m y 
measured at T = 0.4 J are plotted as a function of time for the two cases 1 
(upper panel) and 3 (lower panel). 



we can expect that magnetization at Tsrt heavily fluctuates 
with time. 

To verify weather this reasoning is correct or not, we per- 
formed a standard MC simulation at Tsrt and measured the 
amplitude of magnetization m = {m 2 ± + m 2 ) 1 ^ 2 and the z- 
component of magnetization m z . Figure 4 shows the result. 
The data in the case 1 and those in the case 3 are shown in 
Figs. 4(a) and 4(b), respectively. We see that ra-'s in both the 
cases heavily fluctuate with time, as expected. A similar be- 
havior was observed in the other two components of magne- 
tization m x and m y . However, the amplitude of magnetization 
m is always non-zero because m is non-zero in the low free- 
energy region of Fdifr- Especially, in the case 1, m in the low 
free-energy region is almost constant because Fag has a cir- 
cular structure. This is the reason why fluctuations in m are 
rather small in the case 1 . These results show that, at the SRT 
temperature, spins are aligned in the same direction, and the 
direction changes with time. In other words, the system is in a 
superparamagnetic state at Tsrt. Therefore, the pseudogap in 
magnetization emerges around the SRT temperature because 
of dynamical fluctuations in magnetization which are inherent 
in the superparamagnetic state. 

One may consider that, even if the system is in a in-plane 
magnetized state and ray + 0, the time averages of two mag- 
netizations m x and m y become zero if magnetization freely 
fluctuates in the plane of the film. To examine whether such 



fluctuations occur or not, we measured the time evolutions 
of m x and m y at T — 0.4 J for the two cases 1 and 3. Note 
that the system is in a in-plane magnetized state at this tem- 
perature (see Fig. 1). The result is shown in Fig. 5. We see 
that m x and m y heavily fluctuate with time. This result means 
that the four-fold magnetic anisotropy on square lattice, which 
is induced by magnetic dipolar interactions, is not sufficient 
to prevent in-plane fluctuations of magnetization. To prevent 
these fluctuations, some additional anisotropies such as the 
shape magnetic anisotropy and induced magnetic anisotropy 
which stabilize magnetization to a direction in the plane are 
necessary. 

Lastly, it may be worth pointing out the similarities and 
differences between our model and the normal ferromagnetic 
model which involves only exchange interactions, i.e., our 
model with Cd = C u = 0. A wide spread of low free-energy 
region in the magnetization space and large fluctuations in 
magnetization are also observed in the normal ferromagnetic 
model. However, because exchange interactions are isotropic, 
such behavior is observed at any temperatures in this model. 
On the contrary, this behavior is observed only around Tsrt 
in our model. This is the difference between the two models. 
Because the pseudogap emerges only around Tsrt in exper- 
iments, the normal ferromagnetic model is not proper to de- 
scribe the SRT. 
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5. Conclusions 

We investigated the spin reorientation transition (SRT) in 
ultrathin magnetic films by Monte-Carlo simulations. We con- 
centrated ourselves on the case that the lateral size of the 
film is relatively small and the system has a single-domain 
structure. As a result of free-energy measurement, we found 
that the system is in a superparamagnetic state at Tsrt and 
the disappearance (pseudogap) of magnetization around Tsrt 
emerges due to dynamical fluctuations in magnetization. This 
observation is in contrast to that in large ultrathin magnetic 
films that the pseudogap is caused by a static magnetic struc- 
ture with many complex domains. 
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